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Abstract 

In the present paper the temperature-induced metal-insulator transition 
is studied in a generalized Hubbard model with correlated hopping using re- 
cently obtained expression for energy gap. The dependence of energy gap on 
concentration of doubly occupancy leads to increasing energy gap width with 
increase of temperature. Thus narrow-band system can undergo transition 
from a metallic state to an insulating state with the increase of temperature. 

For some values of intra-atomic Coulomb repulsion U and wq (wq is half- 
bandwidth without taking into account of correlated hopping of electrons) we 
find the values of temperature when narrow-band material undergoes transi- 
tion from a metallic state to an insulating state. We show that at given U /w$ 
metal-insulator transition in model with non-equivalent Hubbard sub-bands 
can occur at smaller temperature than in the Hubbard model. It testifies 
on the fact that taking into account of correlated hopping is important for a 
consideration of metal-insulator transition problem. 

The obtained results are applied to the interpretation of the experimental 
data. 



1 Introduction 

It is known that the electron- hole symmetry is peculiar to the Hubbard model Q. 
One of manifestations of this symmetry is the equivalence of the lower and upper 
Hubbard bands. It is the result of an equality of hopping integrals describing both 
"translational" hoping of holes and doubly occupied sites (doublons) and the pro- 
cesses of their paired creation and destruction. Equality of noted hopping integrals 
is caused by neglecting the matrix elements of electron-electron interaction 



in the general Hamiltonian (the matrix elements (1) describe hopping of electrons 
between % and j lattice sites; 0-function is the Wannier function). 

However, theoretical analysis, on the one hand, and available experimental data, 
on the other hand, point out the fact that the Hubbard model generalization by 
taking into account correlated hopping (1) is of principle necessary - ||. In such 
model hopping integrals describing "translational" hopping of holes and doublons 
are different. These hopping integrals also differ from the hopping integral which is 
connected with the processes of paired creation and destruction of holes and dou- 
blons. In consequence of that the lower and upper Hubbard bands are non-equivalent 
(non-symmetric). In recent years similar models have been studied intensively || - 

0- 

Important puzzle arising in an investigation of these models is metal- 
insulator transition problem which is one of the most essential in narrow-band 
physics |T2[ - WM . In this connection special interest is challenged by the observable 



metal-to-insulator transitions in some narrow-band materials with the increase of 
temperature (see, for example [[14]] - [|18|1). 



On the basis of an approach proposed in the papers || [0| we have studied 
metal-insulator transition in a model of narrow-band material with non-equivalent 
Hubbard subbands (so-called "non-symmetric Hubbard model") at half- filling and 



zero temperature in the paper pOfl . The present paper is devoted to a further study 
of metal-insulator transition in a model with non-equivalent Hubbard subbands, in 
particular, an investigation of temperature-induced metal-to-insulator transition. 



2 Results 

We start from the following natural generalization of the Hubbard model at 
half-filling including correlated hopping (1) 0, |3|, |20| 
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where /i is the chemical potential, af a , (a^) is the creation (destruction) operator of 
an electron of spin a (cr =f, {) on 2-site {a denotes spin projection which is opposite 
to a), n ia = af a a ia is the number operator of electrons of spin a on i-site, U is 
the intra-atomic Coulomb repulsion, t , T 1; T 2 are the integrals of electron hopping 
between nearest neighbors, 

Ti = ^2 J(ikjk), T 2 = J (iiij ) ; 

the primes at the sums in Hamiltonian (2) signify that i ^ j. 

Using a generalized mean-field approximation |3], |19| in Green function method 
we obtain for a paramagnetic state the single-particle energy spectrum as pi| 



where k is the wave vector, Ei(k), (^(k)) are the energies of electrons within 
the lower (upper) Hubbard band, d is the concentration of polar states (holes or 
doublons), t^, t^, t'^ are the Fourier transforms of respective hopping integrals 
t = t + Ti, i = t + 2T 2 , t' = t + T 2 ; t and t are terms describing hopping of 
quasiparticles within the lower and upper Hubbard bands (hopping of holes and 
doublons) respectively, t' describes quasiparticle hopping between hole and doublon 
bands (the processes of paired creation and destruction of holes and doublons). 

Note that single-particle energy spectrum (3) gives the exact atomic and band 
limits: if U = and tk = = t'k = £o(k) (neglecting correlated hopping) then 
£1,2 (k) take the band form, if = tk — t'^ — > then we obtain exact atomic limit. 

With the help of single-particle energy spectrum (3) we find the expression to 
calculate the energy gap width (difference of energies between bottom of the upper 
and top of the lower Hubbard bands): 

AE = -(l-2d)(w + w) + \(Q l + Q 2 ), (5) 



Qi = \j[B(w -w)- Uf + (Adzt') 2 , (6) 
Q 2 = ^[B(w-w) + U] 2 + (4dzt>) 2 , (7) 

where w and w are halfwidths of the lower (hole) and upper (doublon) Hubbard 
bands: w = z\t\, w = z\i\ (z is the number of nearest neighbors to a site). 

At given w, w, t', U and change of d (which can be caused by external influences, 
in particular, by temperature) energy gap (5) vanishes when the condition d < do is 
satisfied where do is the root of the equation 

-(1 - 2d)(w + w) + i(Qi + Q 2 ) = 0. (8) 

Then a metallic state (AE < 0) is realized at d < do, an insulating state (AE > 0) 
is realized at d > do- 

If t = t = t' — to (the Hubbard model) then energy gap (5) takes the following 



form [19]: 



AE = -2w(l - 2d) + ^U 2 + (Adw) 2 , (9) 
and vanishes when the condition d < do is satisfied where 

4 = 1 - { f w)2 (2» > U). (10) 

From Eq. (8) we obtain the dependence of do on U/w ratio (Fig. 1). The param- 
eters Ti = 7i/|to|, T 2 = T 2 /\t \ characterize value of correlated hopping. From Fig. 1 
one can see that value of d depends on the parameters of correlated hopping t%, t 2 
(thus on w/w) weakly when U/w is near zero. But with the increase of U/w the 
value of d begins to depend strongly on the parameters Ti, t 2 . Fig. 1 points out 
the extension of region in which narrow-band material is a metal with the increase 
of w/w ratio (at given U/w). 

The concentration of polar states (obtained with the help of the Green function 
((aiania\a1j a ))) is [|Cj 

1 V T n n 1 



with 
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and chemical potential of narrow-band model with non-equivalent Hubbard sub- 
bands is given by the equation 
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de = 0, 



(14) 



Ei^(e), F £ are obtained from respective formulae (3), (4) for E lf2 (k), F^ by sub- 
stitution of tk — > e, tk — ► (t/t)e, t' k — > (t'/t)e. Here we have used the rectangular 
density of states 



-£^-*(k)) = -^ 2 -£ 2 ), 



2w 



(15) 



where #(:r) = 1 if rr > 0, =0 otherwise; N is the number of lattice sites. 



3 Discussion 



The temperature dependence of chemical potential of a narrow-band model with 
non-equivalent Hubbard subbands obtained from Eq. (12) is plotted in Fig. 2. One 
can see that in the considered model in the region of low and normal temperatures 
chemical potential is essentially dependent not only on the parameters w and w 
and also on temperature (in contrast to the Hubbard model where fj, = U/2), and 
what is more with the decrease of temperature chemical potential rapidly increases 
depending on the paramaters of non-equivalence of Hubbard bands t±, r 2 . In high 
temperature region in the proposed model chemical potential tends to U/2 with the 
increase of temperature; really, at T — > oo the probabilities of an electron finding 
within the lower and upper Hubbard bands (independently of their bandwidths 
ratio) are equal. 

At given U, w, w, t' (constant exterior pressure) concentration of polar 
states (11) increases with the increase of temperature. It leads to the fact that sys- 
tem can undergo transition from the state with AE < to the state with AE > 0, 
i.e. metal-to-insulator transition can occur. In this case the results obtained in the 
Hubbard model and those obtained in non-symmetric Hubbard model can be essen- 
tially different (Fig. 3 illustrates it). Let us take for example U/w = 0.9. One can 
see that at T = the energy gap width in both models is AE < (a metallic state). 
With the increase of temperature metal-to-insulator transition does not occur in the 
Hubbard model, in non-symmetric model the values of parameters Ti, r 2 exist at 
which metal-to-insulator transition occurs. 

In case metal-to-insulator transition occurs in both models from Fig. 3 one can 



model. So, for example, for wq = z\to\ ~ 1.05 eV (such bandwidth of N1S2 was 
estimated in paper [16]) in considered model in a paramagnetic state metal-to- 
insulator transition occurs at T pa 280 K for U/wq = 1.94 and T\ = r 2 = 0.01 
(observable the transition temperature of N1S2 is T ~ 280 K at p ~ 3 MPa [0). 
For the same value of U/wq metal-to-insulator transition occurs at T ^940 K when 
r i — T 2 — (neglecting of correlated hopping, the Hubbard model) and at T = K 
when T\ = r 2 = 0.015. If U/wq = 1.98 then transition from a metallic state to 
an insulating state is realized at T pa 290 K for T\ = r 2 = 0; at T = K when 
T\ = T2 = 0.005. Note that at U ~ 2w the temperatures of metal-to-insulator 
transiton found in both models are essentially different; with a deviation from this 
ratio the difference decreases. 

The obtained temperature dependence of energy gap can explain observable tran- 
sition from the state of a paramagnetic metal to the paramagnetic Mott-Hubbard 
insulator state in the (Vi- x Gr x ) 2 3 compound [T3J, |TJJ in NiS 2 [O and in the 



NiS 2 - x Se :c system [|17|, [L8| with the increase of temperature. 

In summary, metal-to-insulator transition observable in the materials with nar- 
row energy bands can be explained on the basis of the proposed approach at realistic 
values of the parameters characterizing non-equivalence of the Hubbard subbands. 
In this way series of results obtained from such consideration are essentially distinct 
from those obtained without taking into account this non-equivalence (the Hubbard 
model). In particular, the temperature of metal-to-insulator transition is essentially 
smaller than in the Hubbard model (this fact agrees to the observable transition tem- 
perature), and chemical potential disctinguishes from the value p = U/2 (obtained 
in the Hubbard model) and is temperature dependent. 
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Figure 1: The dependence of do on U/w ratio: the upper curve corresponds to 
T\ = T2 = (the Hubbard model); the middle curve - t\ = r 2 = 0.1; the lower curve 

- 7~i = r 2 = 0.2. To the left of respective curve there is a metallic phase, to the right 

- an insulating phase. 

Figure 2: The temperature dependence {0 = ksT) of chemical potential \i for 
U/2w = 1: the upper curve corresponds to T\ = r 2 = 0.3; the lower curve - 
Ti — t 2 = 0.1; the straight line corresponds to values of chemical potential in the 
Hubbard model (ti = r 2 = 0). 

Figure 3: The dependence of energy gap on temperature at U/w = 0.9. The 



